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Abstract 

The paper provides a new framework for the description of linearized adiabatic 
lagrangian perturbations and stability of differentially rotating newtonian stars. In 
doing so it overcomes problems in a previous framework by Dyson and Schutz and 
provides the basis of a rigorous analysis of the stability of such stars. For this the 
governing equation of the oscillations is written as a first order system in time. 
From that system the generator of time evolution is read off and a Hilbert space 
is given where it generates a strongly continuous group. As a consequence the 
governing equation has a well-posed initial value problem. The spectrum of the 
generator relevant for stability considerations is shown to be equal to the spectrum 
of an operator polynomial whose coefficients can be read off from the governing 
equation. Finally, we give for the first time sufficient criteria for stability in the form 
of inequalities for the coefficients of the polynomial. These show that a negative 
canonical energy of the star does not necessarily indicate instability. It is still unclear 
whether these criteria are strong enough to prove stability for realistic stars. 



1 Introduction 



The study of oscillations of stars is an important and exciting field of current astrophysics. 
For instance through period-luminosity and period-radius relationships variable stars pro- 
vide important 'yardsticks' for distance measurements in the universe. Their observation 
yield important information about the interior of stars, like the equation of state of the 
matter, which is otherwise hard to obtain. Further, Neutron star pulsations may be a 
source of gravitational radiation detectable for experiments like LIGO, VIRGO and GEO 
600 in the near future. 



On the other hand, it is probably fair to say that there has not been very much work on 
the mathematical foundations of the theory of stellar oscillations. In the non relativistic 
limit there is a well-known framework for the description of oscillations of nonrotating 
stars |J, (see B for a rigorous version). It is important to note that even on that level 
it turns out that the governing operator of the spheroidal oscillations belongs to a class 
of operators which were apparently (apart from special cases considered in [O]) previ- 
ously unconsidered in operator theory J|] . Somewhat surprisingly it also turned out that 



differently to radial oscillations these operators don't have a compact resolvent. Hence 
the well developed perturbation theory for such operators cannot be applied and a corre- 
sponding theory for the new type of operators still has to be developed. For rotating stars 
there is little known about the relevant operators apart from abstract properties (like the 
symmetry, semiboundedness and continuity of associated operators J2^], |22[]), indication 
of a continuous part in the spectrum [^4, |6| and instabilities caused by so called 'r-modes' 



(also called 'quasi-toroidal modes') JT|, 25]. To my knowledge there is no consideration of 
these operators in sufficient detail. Indeed, a large part of the present paper considers the 
more modest first step of such an investigation, namely to identify and determine which 
operators should be considered in that case. 

The governing equation for linearized adiabatic oscillations of a stationary differentially 



rotating perfect-fuid star in an inertial frame (t, x) is [27 
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£ is the Lagrangian displacement vector field, for j e {1,2,3} the symbol ej denotes the 
canonical unit vector in the direction of Xj, Q is the (bounded open) volume of the star, 
and v,p, p, Ti are the velocity field, pressure, density and the adiabatic index functions of 
the background star satisfying the equations of momentum and mass conservation 

(v ■ V) v = - Q Vp + , V ■ {pv) = . (5) 

and an equation of state. In addition to ([IJ) the variation Sp of the pressure has to vanish 
at the boundary dfl of the star, i.e., 

lim(j?r 1 V-0(2/) = (6) 

for all x e dfl. 

The remarkable paper [12] of Dyson and Schutz provides a framework for deciding 



the stability of the solutions of (|]). In the following this paper is referred to as DS. 
Compared to previous frameworks given in [57], P2J the main step forward in that paper 
is the fact that it relates the stability of the system directly to the growth properties of 



the perturbations in time as is usual for nonrotating stars (see, e.g., [0, ||, 0, 0, |25| 



|36fl). Moreover the paper shows that these growth properties are governed by spectral 



properties of the generator of time evolution. This greatly simplifies the stability discus- 
sion. Unfortunately, the approach still has some drawbacks, and in the present paper 
will be given a varied framework which overcomes those problems. Moreover here are 
given for the first time sufficient criteria for stability in the form of inequalities which 
have to be satisfied by the coefficients of an operator polynomial. These criteria show 
that a negative canonical energy does not necessarily indicate an instability of the star. It 
is still unclear whether these criteria are strong enough to prove stability for realistic stars. 

A rough discussion of the approach of Dyson and Schutz is now given. The paper considers 
axisymmetric solutions of the form 

f (t, x) = expiirmp) f m (t, r, 9) , (7) 

where r, 9, tp are spherical coordinates and m G Z. Inserting this ansatz into ([[]) leads 
to an equation of the same structure with induced operators B' m C' m . The index m is 
supressed in the following discussion. A Hilbert space H' (here X) for the data is chosen 
such that, both, B' becomes continuous and antisymmetric and C becomes symmetric. 
In the nonrotating limit H' goes over into the usual Hilbert space used in the stability 
discussion for spherically symmetric stars. A physically reasonable condition on the back- 
ground model is given which leads to a lower bounded C . Assuming that condition C is 
substituted by its so called Friedrichs extension. This is an abstractly defined self-adjoint 
extension which exists for every densely defined linear symmetric and semibounded op- 
erator in Hilbert space (see e.g. [SIJ Vol. II). [] In the standard way the resulting wave 



equation is written as a first order system in time for £(r, #,£) and (d£/dt)(r,0,t). The 
initial value problem of the system is studied. The Hilbert space of the data is chosen as 
H' 2 with the induced 'euclidean' scalar product. However it is noticed that this is physi- 
cally not meaningful, because the scalar product has no physical interpretation and is not 
even dimensionally correct. From the first order system the linear operator T generating 
time evolution is read of and it is shown that its spectrum is equal to the spectrum of a 
quadratic operator polynomial generated by B' and C |33|. Moreover the resolvent 
of T can be given in terms of the inverses of the operator polynomial. That information 
along with estimates on the resolvent of the operator polynomial are used to give an es- 
timate on the spectrum of T. In general that estimate is not strong enough to decide the 
question of stability of the system. From these estimates it is further shown that there 
is a solution of the initial value problem for the system corresponding to elements of the 
domain of T. The uniqueness of the solution is not shown. The authors remark that 
they could not show that T is the generator of a strongly continuous semigroup and as 
a consequence the results of standard semigroup theory could not be used. Finally, the 
completeness of normal modes of the system is discussed. 

From the description the reader might have noticed that in the derivation of these 
results only abstract properties like 'continuity', 'symmetry', 'semiboundedness' and 'self- 
adjointness' of B' and C play a role. This is indeed true and is the reason why that 
approach is called here a 'framework'. The same also applies to the approach here. As a 



1 The importance of choosing a self-adjoint extension of C can be seen in the limit of no rotation. In 
, H it is shown that for polytropic stars with a polytropic index n < 1 there is an infinite number of 
different self-adjoint extensions which all lead to a well-posed initial value problem for the wave equation. 



consequence these frameworks can be used to describe a lot more physical systems than 
stellar oscillations. The main ingredient for such an application is a system of wave equa- 
tions which is second order in time (with or without first order time derivatives) and which 
is not explicitly time dependent. The 'coefficients' in that system can be (not necessar- 
ily local) linear operators with certain abstract properties. For this reason we abandon 
in Section 2 any reference to rotating stars and just consider abstract wave equations of 
type ([I]). Having this in mind might also provide a better understanding of some of the 
statements below. 



After this digression the discussion of DS is continued. The main problem of the 
approach comes from the chosen Hilbert space along with a scalar product which is not 
related to any physical quantity and not dimensionally correct. Of course the latter could 
be remedied by first introducing a dimensionless time coordinate. But experience tells 
that this should not be essential at such an early stage. Also it is known that the use of 
a suitable Hilbert space decides whether semigroup theory can be applied or not. So it 
is very likely that the use of if' 2 is responsible for the fact that semigroup theory could 
not be applied. Indeed a different choice of the Hilbert space will turn out to be the key 
to the results of this paper. 

Another point which was not addressed in DS is the fact that in addition to (U) the 
boundary condition (||) has to be satisfied that the Lagrangian variation d~p of the pressure 
vanishes at the surface of the star. |25], |9| Indeed it has been shown in || ^§ for the limit 
of no rotation that for a polytropic equation of state with polytropic index n < 1 there 
is a infinite number of different self-adjoint extensions of C, which all lead to different 
initial value formulations for the wave equation. Moreover it has been shown that the 
condition of a vanishing 5p at the surface of the star picks exactly one of these self-adjoint 
extensions. Of course the choice of the Friedrichs extension of C, is equivalent to posing 
a boundary condition. But because of the abstractness of this extension it is not obvious 
and has to be investigated whether it is compatible with (|6]). To my knowledge this has 
been shown only for the case of radial oscillations of spherically symmetric stars in || . 
This point will not be pursued any further in this paper. 



The approach in this paper is similar to that of Dyson and Schutz. The point of 
departure is in the choice of the Hilbert space for the initial data of the first order system. 
Here a space Y is chosen , which is in general a proper subspace of H' 2 . Moreover a 
different and dimensionally correct scalar product is chosen. The square of the induced 
norm of the initial data is a positive definite part of the corresponding canonical energy of 
the system. [16, 17] For this C is split into sum of a strictly positive self-adjoint operator 
A and a 'rest' C . Of course such a decomposition is not unique but it can be shown (see 
Lemmas 14 and 15 in Section 2) that trivial rescalings all lead to the same set Y along 
with equivalent norms onF. In particular such changes lead to theories which are related 
by a similarity transformation and hence the outcome of the stability discussion is not 
affected. In general the canonical energy cannot be used as a norm for Y because it is not 
always positive definite. In situations where it is C can be chosen as zero. In the limit of 
no rotation where B' = and the operator C is semibounded the approach here reduces 
to the approach in [|30| Vol. II (see the proposition at the beginning of chapter X.13) for 
classical wave equations. 



A major consequence of the change is that it allows the use of semigroup theory which 
is a standard and well developed tool in particular in the theory of partial differential 
equations (see e.g. [p^,[|Ilfl,PD||, [|19|] , |29|] , []32| and the cited references therein). This 



simplifies the reasoning a lot, because it can be and will be built on those results. In 
particular here the operator G + which corresponds to T in DS generates a strongly con- 
tinuous group of bounded transformations and hence the well-posedness of the initial value 
problem for the first order system follow from abstract semigroup theory. At the same 
time a considerable generalization is achieved. The operator B' (here denoted by iB) can 
be unbounded and not antisymmetric. Moreover C has not to be assumed symmetric. 
The restrictions imposed on these operators are the following. The operator C has to 
be of the form A + C where A is some densely defined and strictly positive self-adjoint 
operator in X and C is a relatively bounded perturbation of the positive square root 
A 1 / 2 of A. In addition B' has to be a relatively bounded perturbation of A 1 / 2 with rel- 
ative bound smaller than 1. Finally, B' has to be antisymmetric or continuous, but not 
necessarily both. All these conditions are trivially satisfied for the case of axisymmetric 
solutions of ([]]) considered by Dyson and Schutz. Whether this generalization is sufficient 
to provide a framework for ([!]) and not only for its axisymmetric form is not yet clear. 
For this it seems necessary that C given by (^) is semibounded and this is still open. 
The reason for considering also more general situations with nonantisymmetric B' and 
nonsymmetric C is that the framework here will also be used in a future paper in the 
stability discussion of the Teukolsky equations on a Kerr background where this is the 
case. |35| A further important advantage of the approach here is that it can be shown (see 
Theorem 3) that the dominating part of G + (but in general not G + itself) is self-adjoint. 
Using perturbation theory this gives important information on the spectrum of G + and 
is also the basis of the proof that G + is the generator of a strongly continuous group 
of bounded transformations. (See Theorem 7) Further it is the basis for another result 
(see Corollary 12) having no counterpart in DS namely the conservation of the 'canonical 
energy' E. On the other hand it turns out that the spectrum of G + is the same as of 
T. In particular that spectrum is given by the spectrum of the same operator polynomial 
C - XB ' + A 2 , A e C (See Theorem) 



A plausible definition for the stability of a rotating star is the following. The system 
is stable if and only if the semigroup T + (t),t e [0, oo) generated by G + is bounded. Note 
that this definition is invariant to similarity transformations and hence not so sensitive 
to changes of Hilbert space like one only invariant under unitary transformations.^ From 
semigroup theory one has then the following. 

1. The system is unstable if G + has a spectral value with real part smaller than zero. 

2. For a stable system the corresponding spectrum of G + is contained in the closed 
right half plane of the complex plane. 

3. From only the fact that the spectrum of G + is part of the closed right half plane of 
the complex plane, it cannot be concluded that the system is stable. Q 

4. The system is stable if the real part of all 'expectation values' 

m + o (8) 

2 Such a definition would be given for instance by the demand that the semigroup should be contractive, 
i.e., that the norms of the semigroup elements are ^ 1. 

3 For a counterexample compare for instance the note after the proof of Corollary 9. 



is positive 0) for all elements £ from the domain (or a core) of G+.f] 

Point 1 gives a sufficient but not necessary condition for instability. Note that this con- 
dition is invariant under similarity transformations. Moreover because of Theorem 13 it 
is equivalent to the condition that there is complex number A with real part smaller than 
zero such that 

C'-XB' + X 2 (9) 

is not bijective. This reduces in the nonrotating case to the condition that C is strictly 
negative, which is a well known sufficient condition for instability.pl An important final 
observation is that from the existence of such a A follows the existence of an element 
£ from the Hilbert space such that the corresponding function of norms \T + (t)£\,t G R 
grows exponentially for large times.Q Hence the existence of such a A leads to a much 
stronger kind of instability. 

Point 4 gives a sufficient but not necessary condition for stability. It is appealing be- 
cause it is of the form of an inequality, which is more easily accessible than the spectrum 
of G + . On the other hand it is strong and not invariant to similarity transformations. 
It turns out to be equivalent to C being strictly positive, i.e., that the spectrum of this 
operator consists only of positive real numbers different from zero. Note that this reduces 
to a known sufficient condition for stability in the nonrotating case. But for such stars 
it can be satisfied only for radial oscillations (for instance this is the case for constant 
Ti > 4/3), but not for nonradial oscillations.^ |2|, 0] Note that in the limit of no rotation 
the trivial toroidal oscillations give rise to solutions of (|l]J|) whose norm increases linear 
in time for large times. Hence applying the stability definition above to that limit would 
lead to an 'unstable star'. Of course, these oscillations can be excluded in that case just 
by considering a reduced operator. 

The following two stability criteria are new. They will be proven in Theorem 17. 

1. If B' and C are such that 

<i\C'i>-- A <i\B'i> 2 ^Q (10) 

for all £ from the domain of C such that ||£|| = 1 then the spectrum of G + is purely 
imaginary. 

2. If the operator 

ib b 2 

C'-™B'-- (11) 
2 4 K J 

is positive for some b G R then the spectrum of G + is purely imaginary and there 

are K ^ and t ^ such that 



\T + (t)\^Kt (12) 



for all t > t 



o- 



4 Then G+ generates a contraction semigroup. 

5 This is easily seen for instance by using Theorem 4.1 in chapter 4 of ^9|. Here it is important 
to remember that in general the spectrum of G+ does not only consist of 'eigenvalues' (for which this 
statement is of course trivially satisfied) but also values \i S C for which the map G+ — /i is just not onto. 
Such values are often from a continuous part of the spectrum. 

6 This is obvious since the spectrum of the trivial toroidal oscillations is {0}. But in || it has also 
been shown that is in the spectrum of spheroidal oscillations. 



Note for the first point that —(1/4) < £\B'£ > 2 is positive, because of the antisymmetry 
of B' . Also note in this connection that in DS it has been shown that C — (l/4)5 /2 
is bounded from below uniformly in m. Unfortunately, in general this does not imply (|i~0p. 

It is still unclear whether these criteria are strong enough to prove stability for realistic 
stars. On the other hand the second criterium has been sucessfully applied in the stability 
discussion of the Kerr metric where the master equation governing perturbations is of the 
form (D, too.§ 



2 The framework 

This section developes the initial value formulation for abstract differential equations of 
the form (|l|). It is self-contained and necessarily very technical. The reader who is not 
interested in the excessive mathematical details given here is referred to the introduction. 



The used nomenclature can be found in standard textbooks on Functional analysis. [30] 
Vol. I, ill M 



Before going into the mathematical details it is explained about the meaning of the 
individual results of this section. The section is based on the assumptions General As- 
sumption 1 and General Assumption 4 on three operators A, B and C. A different form 
of General Assumption 1 which is more convenient for applications can be given in the 
obvious way using Lemma 18. Definition 2 gives the Hilbert space Y which is used here 
instead of the Hilbert space in DS. A rigorous form ( |53D of ([]]) along with the existence 
and uniqueness of the solution corresponding to initial values is given in Theorem 11. 
Corollary 12 gives the corresponding 'energy' along with an identity for its time deriva- 
tive. The analogue G + of the generator T in DS is given in Definition 5. Theorem 3 proves 
that the 'dominating parts' of G + are self-adjoint. In Theorem 7 it is proved that under 
General Assumption 1 and General Assumption 4, both, G + and — G + are generators of 
strongly continuous semigroups T + and T_, resp. Theorem 13 shows the identity of the 
spectrum of G + with the spectrum of an operator polynomial generated by the operators 
B and A + C . Lemmas 14 and 15 show that certain simple rescalings of A and C which 
formally leave invariant lead to theories which are related by a similarity transfor- 
mation. Theorem 16 shows for a special case how these rescalings can be used to derive 
a better estimate for the growth of T + and T_ than the one induced by ([$(]) in Lemma 6. 
Theorem 17 gives sufficient criteria for stability in the form of inequalities which have to 
be satisfied by the coefficients of the operator polynomial. Part (ii) of this Theorem has 
been sucessfully applied in the discussion of the stability of the Kerr metric. || 

The rest of this section contains the mathematical details. 

Assumption 1 In the following let (X, < | >) be a non trivial complex Hilbert space. 
Denote by \\ \\ the norm induced on X by < \ >. Further let A : D(A) X be a densely 
defined linear self-adjoint operator in X for which there is an e G (0, oo) such that 

<^\AC>^e<^> (13) 

for all £ G D(A). Denote by A 1 ' 2 the square root of A with domain D(A l l 2 ). Further let 
be B : D(A 1 ^ 2 ) — > X a linear operator in X such that for some a 6 [0, 1) and & 6 I 

\\Ba 2 ^a 2 \\AV 2 a 2 + b 2 U\\ 2 (14) 



for all £ G D(A 1 / 2 ). Finally, let C : D{A 1 / 2 ) — > X be linear and such that for some real 
numbers c and d 

\\ca 2 ^ c'wa^w 2 + d 2 u\\ 2 as) 

for all£e D(A l l 2 ). 

Note that as a consequence of (|l^) the spectrum of A is contained in the interval [e, oo). 
Hence A is in particular positive and bijective and there is a uniquely defined linear and 
positive selfadjoint operator A 1 / 2 : D(A l l 2 ) — > X such that (A 1 / 2 ) 2 = A. That operator is 
the so called square root of A. Further note that from its definition and the bijectivity of A 
follows that A 1 / 2 is in particular bijective. This can be concluded for instance as follows. 
By using the fact that A x l 2 commutes with A it easy to see that for every A G [0, e 1 ^ 2 ) by 
{A 1 / 2 + \)(A — A 2 ) -1 there is given the inverse to A 1 ! 2 — A. Hence the spectrum of A 1 ! 2 
is contained in the interval [e 1 ^ 2 , oo). All these facts will be used later on. 

Definition 2 We define 

Y := D(A 1 ^ 2 ) x X (16) 

and ( | ) : Y 2 -> C by 

:=< A^tlA^m > + < 61*72 > (17) 
for all £ = £ 2 ), 77 = (771,772) G K. 
Then we have the following 

Theorem 3 (i) (Y, ( | )) is a complex Hilbert space. 

(ii) The operator H : D(A) x D(A 1/2 ) —>Y inY defined by 

F£:=H6,^6) (18) 
for all £ = (£1,6) e -D(^) x -D(^4 1//2 ) is densely- defined, linear and self-adjoint. 

(iii) The operator B : D(H) — > Y defined by 

Bi := (0, -56) (19) 

for all £ = (656) £ D(H) is linear. If B is symmetric then B is symmetric, too. 
If B is bounded then B is bounded, too, and the corresponding operator norms \\B\\ 
and \B\ satisfy 

\B\ ^ ||B|| . (20) 

(iv) The sum H + B is closed. If B is symmetric then H + B is self-adjoint. 

(v) The operator V : Y — > Y defined by 

V£:=(0,iC£i) (21) 
for all £ = (6)6) G V is linear and bounded. The operator norm \V\ ofV satisfies 

\V\^{c 2 + d 2 /ef' 2 . (22) 



Proof: (i): Obviously, ( | ) defines a hermitean sesquilinear form on Y 2 . That ( | ) is 
further positive definite follows from the positive definiteness of < | > and the injectivity 
of A 1 / 2 . Finally, the completeness of (Y, | |), where | | denotes the norm on Y induced 
by ( | ), follows from the completeness of (X, || ||) together with the fact that A 1 ! 2 has 
a bounded inverse. Here it is essentially used that is not contained in the spectrum 
of A. (ii): That D(A) x D(A 1/2 ) is dense in Y is an obvious consequence of the facts 
that D(A) is a core for A 1 / 2 (see e.g. Theorem 3.24 in chapter V.3 of |2~3|) and that 



D(A 1 / 2 ) is dense in X. The linearity of H is obvious. Also the symmetry of H follows 
straighforwardly from the symmetry of A 1 / 2 . By that symmetry one gets further for any 
i — (Ci)^) G D(H*) and any 77 = (771,772) G D(H): 

(H*£\ V ) = <(H*0 1 \Ar h > + <(H*0 2 \V2> 
= ^\HV) = <-t&\A Vl > + <tA 1 / 2 Z 1 \A 1/2 V2> (23) 

and from this by using that A is bijective and A 1 / 2 is self-adjoint that £1 G D(A) and 

(iTOi = -*6 , (H*0 2 = %Mi ■ (24) 

Hence H is an extension of H* and thus H = H*. (iii): The linearity of B is obvious. 
Also it is straightforward to see that B is symmetric if B is symmetric. If B is bounded 
then 

\m 2 = \\B&\\ 2 < ||5|| 2 ||6|| 2 ^l|5|| 2 |£| 2 (25) 

for all £ = (^1,^2) G D(H). Hence B is also bounded and \B\, \\B\\ satisfy the claimed 
inequality, (iv): Obviously, (|TJ]) implies 

\B£\ 2 ^ a 2 \H£\ 2 + b 2 \£\ 2 (26) 

for all £ G D(H). From this it is easily seen that H + B is closed (see, e.g., |Tj|, Lemma 
V.3. 5). Moreover in the case that B (and hence by (iii) also B) is symmetric (p6|) implies 
according to the Kato-Rellich Theorem (see, e. g., Theorem X.12 in |3(| Vol. II) that 



H + B is self-adjoint. For the application of these theorems the assumption a < 1 made 
above is essential, (v) The linearity of V is obvious. For any £ = (£1,^2) G V one has 

\V£\ 2 = i^c^w 2 ^ c 2 \\A^ 2 ^\\ 2 + d 2 ]]^]] 2 

= c 2 \\A l / 2 ^\\ 2 + d 2 \\(A 1 / 2 )- 1 A 1 /%\\ 2 < (c 2 + d 2 / e )|£| 2 . (27) 

In the last step it has been used that 

|| (A 1 / 2 ) -1 1| ^ 1/y/e . (28) 

This follows by an application of the spectral theorem (see, e.g. Theorem VIII. 5 in [0O[ 
Vol. I) to A 1 ! 2 . Since £ is otherwise arbitrary from ( |27| ) follows the boundedness of V 
and the claimed inequality. □ 

Assumption 4 In i/ie following we assume in addition that B is symmetric or bounded. 
Note that condition (|14]) is trivially satisfied if B is bounded. We define: 
Definition 5 

G + := -i(H + B + V) , G_ :=i(H + B + V) . (29) 



then 

Lemma 6 The operators G + and G- are closed and quasi- accretive. In particular 

Re(Z\GO>-{»B + \V\)(t\0 (30) 

for G G {G + , G_} and all £ G D(H). Here Re denotes the real part and 

( if B is symmetric , , 

^ B := \ ||S|| ifB is bounded ' ^ > 

Proof: That G + and G_ are closed is an obvious consequence of (iv) and (v) of the 
previous theorem. Further if B is symmetric one has because of (iv) and (v) of the 
preceeding theorem 

Re(t\G±o = ^Re(c\m) > > -\v\ m m 

for all £ G D(H). Similarly, if B is bounded one has because of (ii),(iii),(iv), (f20|) 

MZ\G±£) = TRe(Z\i(B + V)£) > -\(Z\i(B + V)£)\ > -(\\B\\ + \V\) (£|£) (33) 

for all £ G D(H). Hence in both cases G + and G_ are quasi-accretive. □ 

Theorem 7 The operators G + and G_ are infinitesimal generators of strongly continuous 
semigroups T + : [0, oo) — > L(Y, Y) and T_ : [0, oo) — > L(Y, Y), respectively. If (jl± G R are 
such that 

#e(£|G ± £) > -li± (£|£) (34) 

for all £ G D(H) the spectra of G + and G_ are contained in the half-plane [— oo) x R 
and [— /i_,oo) x R, respectively, and 



\T+(t)\ ^ exp(/x+t) , |T_(*)| ^ exp(M) (35) 



/or a// i G [0, oo). 



Proof: Obviously, by the Lumer-Phillips theorem (see, e.g., Theorem X.48 in Vol. II of 
p0| ) and the preceeding lemma the theorem follows if we can show that there is a real 
number A < min{— //+, — //_} such that G± — A has a dense range in Y . For that proof 
let be £ some element of D(H) and A any real number such that |A| ^ \V\ 2 . Then we get 
from the symmetry of H 

\(H - tX)^ 2 = \Htf + \ 2 \tf (36) 

and 

\(H- i \)Z\>m ax {\HZ\,\\\ 1/2 \VZ\} . (37) 
Using these identities together with ( |14"D 



|(5 + T0£| 2 < \m 2 + 2\BZ\\VZ\ + \VZ\ 2 

< a 2 |tf£| 2 + 2|5£||V£| + (6 2 +|\/| 2 )|£| 2 
^ a 2 |tf£| 2 + 2a|tf£||U£| + (&+|V|) 2 |£| 2 

< a 2 |(# - *A)£| 2 + 2a\(H - zA)£| |F£| + [(& + M) 2 - « 2 A 2 ]|£| 2 

< a(a + 2|A|- 1 / 2 )|(if-zA)£| 2 + [(6+|\/|) 2 -a 2 A 2 ]|£| 2 . (3* 



Hence for any real A with 

|A| > max{|Vf ,4(1 -a)- 2 ,(6+ |\/|)/a, M} , (39) 

where we assume without restriction that a > 0, we get 

\{B + V)t\^a'\{H-i\)t\ (40) 

where a 1 is some real number from [0,1). Since £ G D(H) is otherwise arbitrary, we 
conclude that 

(B + V)(H -iX)- 1 (41) 
defines a bounded linear operator on Y with operator norm smaller than 1. Since 

H + B + V - iX = (l + (S + V){H - iX)-^j (H - iX) (42) 

we conclude that H + B -\-V — iX is bijectice and hence also that G + — X and G_ — A are 
both bijective. Hence the theorem follows. □ 

We note that General Assumption 4 has been used only to conclude that G + and G_ 
are both quasi-accretive. Now it is easy to see that if B is in addition such that iB is 
quasi-accretive (but not necessarily bounded or antisymmetric) then — iB and hence also 
G + are quasi-accretive, too. As a consequence we have the following 

Corollary 8 Instead of General Assumption 4 let B be such that iB is quasi- accretive. 
Then G + is the infinitesimal generator of a strongly continuous semigroup T + : [0, oo) — > 
L(Y,Y). Iffi+eR is such that 

Re(£\G + > -n + m (43) 
for all £ G D{H) the spectrum of G + is contained in the half-plane [—//+, oo) x R and 

\T + (t)\ ^ exp(/M) (44) 

for all t G [0, oo). 

Theorem 7 has the following 

Corollary 9 (i) By 

nt):=\ If I S ° T ?l (45) 
v ' \ T-(-t) fort < v ' 

for allt EM. there is defined a strongly continuous group T : 1R — > L(Y, Y). 

(ii) For every t G K. and every £ G D(G + ) there is a uniquely determined differentiable 
map u : H. — > F sfic/i i/iat 

«(<o) = C (46) 

and 

u'(t) = -G + u(t) (47) 
/or allt G 1R. i/ere ' denotes differentiation of functions assuming values in Y. 



(iii) The function (u\u) : R — > R defined by 

(u\u)(t) := (u(t)\u(t) ) ,t G R (48) 

is differentiable and 

(u\u)'(t) = -2Re ( u(t) ) (49) 

/or all t G R. 

Proof: The corollary follows from Theorem 7 by standard results of semigroup theory. 
For instance, see section 1.6 in [29] for (i) and section IX. 3 in [23] for (ii). (iii) is an 
obvious consequence of (ii). □ 

Note in particular the the special case [] that there is a non trivial element rj in the 
kernel of A + C for which there is £ G D(A) such that 

(A + C)£ = -iBr] . (50) 

Then by 

u(t) := (£ + tr),r)) , t G R (51) 



there is given a growing solution of (4"?D . 

The following lemma is needed in the formulation of the subsequent theorem. 
Lemma 10 By 

HeiLv. := ll^ll , £ G -D(v4 1/2 ) (52) 
i/iere zs defined a norm || H^i/2 on D(A 1 ^ 2 ). Moreover 

Wi := (-D(A 1/2 ), || |U/ 2 ) (53) 

complete. 

Proof: The lemma is a trivial consequence of the completeness of X and the bijectivity 
oiA^.n 

Theorem 11 Let be to G R ; £ G -D(A) and n G -D^ 1 / 2 ). Then there is a uniquely 
determined differentiable map u : R — > W\ with 

u(to) = £ and u'(t ) = rj (54) 

and such that u' : R — > X is differentiable with 

(u')'{t)+iBu'(t) + (A + C)u(t) = (55) 

for all t G R. 

7 Such cases are easy to construct. 



Proof: For this let be v — (vi, v%) : R — > Y be such that 



v(t ) = (£,r)) (56) 

and 

v '(t) = -G+v(t) ,t e R . (57) 

Such u exists according to Corollary 9 (ii). Using the continuity of the canonical projec- 
tions of Y onto W\ and X it is easy to see that u :— V\ is a differentiable map into VFi 
such that u' : R — > X is differentiable and such that (|54|), (|55|) are both satisfied. On 
the other hand if u : R — > W\ has the properties stated in the corollary it follows by the 
continuity of the canonical imbeddings of W%, X into Y that w := (u, u') satisfies both 
equations (|56D and (|57|) . Then u = V\ follows by Corollary 9 (ii). □ 

Corollary 12 In addition to the assumptions made let C be in particular bounded. f\ 
Further let u : R — >• W\ be differentiable with a differentiable derivative u' : R — >• X and 
such that (5t) holds. Finally, define E u : R — » R by 

E u {t) := ]- (< u\t)\u\t) > + < u{t)\{A + Re(C))u(t) >) . (58) 

Then E u is differentiable and 

jpt/.\ _ / —Im < u{t)\Im{C)u' {t) > for symmetric B , . 

= I \ < u'{t)\Im(B)u'{t) > -Im < u(t)\Im{C)u'{t) > for bounded B (iJ ' : ' 

for all t 6 R, where for any bounded linear operator F on X: 

Re(F) :=\{F + F*) , Im(F) := 1 (F - F*) . (60) 

Proof: For this define t> := (u,u f ). Then according to the preceeding proof v satisfies 
(|57|). For a symmetric B it follows by Corollary 9 and Theorem 3 (iv) that 

(v\v)'(t) = 2Re(v(t)\iVv(t)) 

= - < u'(t)\Cu{t) > - < Cu{t)\u'(t) > (61) 
= - < u\Re(C)u>' (t) - 21m < u(t)\Im{C)u\t) > 

for all t G R. In the last step it has been used that u is also differentiable with the 
same derivative viewed as map with values in X. This follows from the fact the canonical 
imbedding of W\ into X is continuous since A 1 / 2 is bijective. Further the definition 

< u\Re(C)u > (t) :=< u{t)\Re{C)u{t) > , t e R (62) 



for the map < u\Re(C)u >: R — ► R has been used. Obviously, ( p9|) follows from ( pT|) by 
using definition ([581) . In this step also the symmetry of A 1 / 2 is used together with the 
fact that u assumes values in D(A). For a bounded B by Corollary 9 and Theorem 3 (ii) 
follows that 

(v\v)'(t)= 2Re(v(t)\i(B + V)v(t)) 
2Im<u'(t)\Bu'(t) > 

- < u'(t)\Cu(t) > - < Cu{t)\u'(t) > (63) 
2 < u'(t)\Im(B)u'(t) > 

- < u\Re(C)u>' (t) - 21m < u(t)\Im(C)u'(t) > 



'Note that in this case (fL5[) is trivially satisfied 



for all t G R. Obviously, fl59|) follows from (51) by using definition 



□ 



The next theorem relates the spectrum of G + to the spectrum of the so called opera- 
tor polynomial A + C — XB — A 2 , where A runs through the complex numbers. [28, j33fl 



Theorem 13 Let A be some complex number. 

(i) Then H + B + V — A zs noi injective if and only if A + C — XB — X 2 is not injective. 

IfH + B + V — X is not injective then 

ker(H + B + V - X) = {(£, zA£) : £ G ker(A + C - Afi - A 2 )} . (64) 

(ii) Further H + B + V — A zs bijective if and only if A + C — AS — A 2 is bijective. If 

H + B + V — X is bijective then for all rj = (rji, r/ 2 ) G Y: 

(H + B + V- A)- 1 // = (f , *(A£ + m)) , (65) 

w h ere 

i = (A + C -XB - X 2 )- 1 ^ + A)f7i - ir} 2 ) . (66) 

Proof: (i) If H + B + V - A is not injective and £ = £ 2 ) G ker(iJ + B + V - X) it 
follows from the definitions in theorem 3 that 

6 = zA£l , (A + C - XB - A 2 )£i = (67) 

and hence also that A + C — XB — X 2 is not injective. If A + C — XB — X 2 is not injective 
it follows again from the definitions in theorem 3 that 

(H + B + V-X)(Z,iX£) = (68) 

and hence also that H + B + V — A is not injective. (ii) If H + B + V — A is bijective 
it follows by (i) that A + C — XB — X 2 is injective. For r] G X and £ = (£1,^2) : = 
(H + B + V — A) _1 (0, iry) it follows from the definitions in theorem 3 that 

(A + C -XB- A 2 )£i = 77 (69) 

and hence that A + C — XB — X 2 is also surjective. If A + C — A5 — A 2 is bijective it 
follows by (i) that H + B + V — A is injective. Further if rj = (rji, 772) G K and £ is defined 
by (|S6|) it follows from the definitions in theorem 3 that 



(ff + 5 + V - X)(t *(A£ + m)) = V (70) 

and hence that H + B + V — A is also surjective. □ 
Lemma 14 Le£ be e' < e and 

A' := A - e' , C :=C + s' . (71) 



Then 



(i) 

D(A n / 2 ) = D^A 1 ' 2 ) (72) 

and for all £ G _D(v4 1/2 ) 

\\A 1 / 2 £f = WA'^if + e'Uf . (73) 

(ii) The operators A',B and C satisfy 

<£\A'Z> > ( £ _ £ ')<£|£> 

\\BZ\\ 2 < a 2 \\A^ 2 i\\ 2 + {a 2 e' + b 2 )\\i\\ 2 
llC'eil 2 < |c| [\c\ + 2\e'\{e-e')- 112 ] \\A' l '\f + 

[|e'| + (c 2 |e'|+^ 2 ) 1/2 ] 2 lief (74) 

/or a/Z£ G £>(A 1/2 ). 

Proof: (i) First, since e' < e by (|71~|) there is defined a linear self-adjoint and positive 
operator A' in X. Obviously, using the symmetry of A x l 2 and A' 1 / 2 (|73[ ) follows for all 
elements of D(A). From this (|72|) and (|73| ) follow straightforwardly by using the facts 
that D(A) is a core for both, v4 1//2 and A' x l 2 (see e.g. Theorem 3.24 in chapter V.3 of [p3| ). 
that X is complete and that both operators, A 1//2 and A' 1 / 2 are closed, (ii) The first two 
inequalities are obvious consequences of the corresponding ones in General Assumption 
1, the definition ([n]) and of (|73|). For the proof of the third we notice that from the first 
inequality along with an application of the spectral theorem (see, e.g. Theorem VIII. 5 in 
|3T| Vol. I) to A' 1/2 follows that 

||(A /1 / 2 )- 1 Kl/x/^ 7 . (75) 

Further from General Assumption 1 and ([F^) one gets 

\\Ca 2 ^c 2 \\A n / 2 a 2 + (c 2 \e'\+d 2 )U\\ 2 ■ (76) 
for all £ G D(A l l 2 ). From these inequalities we get 

HC'eir < \\Ci\\ 2 + 2\e'\\\Ci\\U\\+e' 2 U\\ 2 (77) 
*C c 2 \\A n / 2 a 2 + (e' 2 + c 2 \e'\ + d 2 )U\\ 2 + 2\e'\ ||C£|| ||f || 
sC c 2 \\A' l l 2 i\\ 2 +[\e'\ + (c 2 \e'\ + d 2 ) 1 ' 2 ] 2 ||£|| 2 + 2\e'\ \c\ \\A n ^\\ 
*C \c\ [\c\+2\e'\{e-e'r l/2 ] U' l ' 2 i\\ 2 + [ \e'\ + (c 2 \e'\ + d 2 ) 1 / 2 ] 2 ||£|| 2 
for all £ G D(A l l 2 ) and hence the third inequality. □ 

As a consequence of (ii) the sequence X, A',B,C satisfies General Assumption 1. The 
corresponding Y given by Definition 2 is because of (i) again given by (0). Moreover the 
corresponding norm | | ' on Y turns out to be equivalent to | | . More precisely one has for 
every e' ^ 

Lemma 15 

\\'^ \ l^e^is-s')-^ 2 ] I' (78) 
and for every bounded linear operator F on Y : 

e~ 1 ' 2 (e - e 1 ) 1 ' 2 \F\ ' < |F| < e 1 ' 2 (e - e')- 1 ' 2 \F\ ' . (79) 



Proof: The first inequality is a straightforward consequence of and ([75]). The second 
inequality is a straightforward implication of the first. □. 

Note that the G± corresponding to the the sequence X, A', B,C are the same for all 
e' (e' drops out of the definition). Moreover as a consequence of the preceding lemma the 
the topologies induced on Y are equivalent. Hence the generated groups are the same, 
too. This will be used in the following important special case. 

Theorem 16 Let be A = A + e, where A is a densely defined linear positive self-adjoint 
operator and let be C = —e. Then 

\T±(t)\ ^ ee 1/2 t exp{fi B t) (80) 

for allt^ e- 1 ' 2 . 

Proof: For this let be e' G [0,e) and define A' and C as in Lemma 14. Hence 

A' = A + e-e' , C = -{e - e') . (81) 
Then from Theorem 3(v), Lemma 6 and Theorem 7 we conclude that 

|T ± (i)| / ^ex P ([/i B + ( £ -e0 1/2 ]O (82) 
for all i 6 R and hence by Lemma 15 that 

|T ± (t)|^e 1 / 2 ( £ - £ ')" 1/2 exp([ /UB + (5- £ ') 1/2 ]0 • (83) 
For t ^ e~ l l 2 we get from this fl80D by choosing 

e':=e-t- 2 . u (84) 
Note that in this special case (^) is conserved and positive. 
We are now giving stability criteria. 

Theorem 17 In addition let B and C be both symmetric. 
(i) Let A, B and C be such that 

<^A + C^>+^<^> 2 ^0 (85) 

for all £ G D(A) with ||£|| = 1. Then the spectrum ofiG + is real. 

(ii) In addition let B and C be both bounded and let A + C + (b/2)B — (6 2 /4) be positive 
for some b G R. Then the spectrum of iG + is real and there are K ^ and to ^ 
such that 

\T(t)\^K\t\ (86) 

for all \t\ ^ to- 



Proof: (i): First from General Assumption 1 and the assumed symmetry of B and C 
follows that, both, by A~ X I 2 BA~ X I 2 and A~ l l 2 C A~ 1 / 2 there is given a bounded symmetric 
and hence (by the theorem of Hellinger and Toplitz) also self-adjoint operator on X . 
Hence 

A(X) := X 2 A~ 1 + XA- 1/2 BA- 1/2 - (l + A~ 1/2 CA^ 2 ) , A G C (87) 

defines a self-adjoint operator polynomial in L(X,X). In addition one has A^ 1 ^ 1/e . 
Further for every f G D(A 1/2 ) and A G C 

< £|A(A)£ >=< rj\A 1/2 A(\)A 1/2 r) >= - < r)\(A + C - XB - X 2 )t] > (88) 

where r\ := A~ x l 2 ^ G D(A). Now (|B5|) implies that the roots of the polynomial < r)\(A + 
C - XB - X 2 )r] >, A G C are real. Hence by (|||) the roots of < £|A(A)£ >, A G C are 
real, too. Since £ G D(A 1 ^ 2 ) is otherwise arbitrary and D(A 1 ^ 2 ) is dense in X this implies 
also that < £|A(A)£ > has only real roots for all £ G X. Hence (see p8| , Lemma 31.1) 



the polynomial A(X), X G C is weakly hyperbolic and has therefore a real spectrum. As a 
consequence A(X) is bijective for all non real A. Now for any such A 

A + C-XB-X 2 = -A 1 ' 2 A i (A)A 1 / 2 , (89) 

where A 1 / 2 denotes the restriction of A 1 / 2 , both, to D(A) in domain and D(A 1 ^ 2 ) in range 
and A (A) denotes the restriction of ^4(A) to D(A 1 ^ 2 ), both, in domain and in range. For 
this note that A(X) leaves D(A 1 ^ 2 ) invariant. Further from the bijectivity of A 1 ^ 2 , A(X) 
and ( |37"D follows the bijectivity of A 1 / 2 and A _ (A), respectively and hence by (|89|) that 
A + C — XB — X 2 is bijective. This is true for all non real A and hence it follows by 
Theorem 13 that the spectrum of iG + is real, (ii) So let B and C be both bounded and 
let A + C + (b/2)B — (b 2 /4) be positive for some b G R. In addition let be e some real 
number greater than zero and define 

A' := A + C + (b/2)B - (6 2 /4) + e , C := -e , B' := B — b . (90) 

First it is observed that 

D(A ,1/2 ) = D(A 1/2 ) (91) 

and that there exist nonvanishing real constants K\ and K 2 such that 

KlWA 1 ' 2 ^ 2 ^ \\A'V 2 ti\\ 2 < tff IIAV^II 2 (92) 

for every £ G /^(A 1 / 2 ). This can be proved as follows. Obviously, by the symmetry of A 1 / 2 
and A' 1 / 2 , the Cauchy-Schwarz inequality, the boundedness of B, C, A~ x l 2 and A'~ x l 2 
follows the existence of nonvanishing real constants K\ and K 2 such that (|9~2"D is valid for 
all £ G .D(.A). Since D(A) is a core for both, A 1 / 2 and A' 1 / 2 (see e.g. Theorem 3.24 in 
chapter V.3 of (23|) from that inequality follows (0) and (0) for all £ G D(A 1/2 ). Note 
that to conclude this it is used that X is complete and that, both, A 1 / 2 and A' x l 2 are 
closed. 

Obviously, from the assumptions made follows that also A',B' and C instead of A, 
B and C, respectively, satisfy General Assumption 1 and General Assumption 4. Hence 
by Theorem 16 follows that 

\Tl(t)\' ^ ee 1/2 t (93) 



for all t ^ e~ 1//2 , where primes indicate quantities whose definition uses one or more of the 
operators A', B' and C instead of A, B and C. In addition (|9l| ) and ( p2|) imply K = Y' 
as well as the equivalence of the norms | | and | | '. Now define the auxiliar transformation 
So : Y> -h. Y by 

So£:=(£i,£ 2 -z(&/2)£i) (94) 

for all £ = (£1,^2) ^ Obviously, So is bijective and bounded with the bounded inverse 
Sq" 1 given by S^f := (£1, £2 + W 2 )£i) for all £ = (£ x ,£ 2 ) G K. In addition define 
S ± :[0,oo)^L(F,F)by 

S ± (t) := exp(Ttbt/2)S Tl(t)S^ , (95) 

for all t G [0, 00). Obviously, S± defines a strongly continuous semigroup with the corre- 
sponding generator 

SoG^So^i^Gi . (96) 



This implies S± = T± and by fl93|) and fl9q) the existence of K ^ and to ^ such that 
flSEp is valid for all \t\ ^ t . Finally, from this follows by the Theorem of Hille-Yosida- 
Phillips that the spectrum of iG + is real.n 

Lemma 18 Let D be a core for A. Further let be Bq : D — > X a linear operator in X 
such that for some real numbers a and 6 

\\B a 2 ^al<Z\AZ>+b 2 U\\ 2 (97) 

for all £ G D. Then there is a uniquely determined linear extension Bo : D(A 1 ^ 2 ) — > X of 
Bq such that 

\\M\ 2 <$U 1/2 t\\ 2 + %\\t\\ 2 ( 98 ) 

for all £ G D(A l l 2 ). If Bq is in addition symmetric B is symmetric, too. 

Proof: First we notice that D is a core for A 1 / 2 , too. Obviously, since D(A) is a core for 
A 1 / 2 (see e.g. Theorem 3.24 in chapter V.3 of f23|) this follows if we can show that the 
closure of the restriction of A 1 / 2 to D extends the restriction of A 1 / 2 to D(A). To prove 
this let £ be some element of D(A). Since D is a core for A there is a sequence £ , £1 • • • 
of elements of D converging to £ and at the same time such that A£ , A£x . . . converges to 
A£. Since A l/2 has a bounded inverse it follows from this that A 1/2 £ , A 1/2 £ 1 . . . vonverges 
to A 1//2 £. Since £ can be chosen otherwise arbitrarily it follows that the closure of the 
restriction of A 1 / 2 to D extends the restriction of A 1 / 2 to D(A) and hence that D is a core 
for A 1 / 2 . Hence for any £ G D^A 1 ^ 2 ) there is a a sequence £o,£i • • • in D converging to £ 
and at the same time such that A 1//2 £o, A 1 / 2 ^ ... is converging to A 1 / 2 ^. Hence by (|97|) 
along with the completeness of X follows the convergence of the sequence -B £i, -Bo£2 • • • 
to some element _B£ of X and 

\\Ba 2 ^a 2 \\A^ 2 a 2 + b 2 oU\\ 2 ■ (99) 



Moreover if £ ', £/ . . . is another sequence having the same properties as £0, £1 • • • by (J9 
follows that 



5£ = lim E £ n = lim B & 



(100) 



From this it easily seen that by defining 



B:={p(A x l 2 )^X,Z»B£) (101) 

there is also given a linear map. Hence the existence of a linear extension of B satisfying 
(El) is shown. Moreover from the definition it is obvious that B is symmetric if Bq is in 



addition symmetric. If on the other hand Bo is a linear extension of B$ satisfying (|98 
and £ and £i,£2 are as above from (p8|) follows that 



B £ = lim B £n ■ (102) 
Finally, since £ can be chosen otherwise arbitrarily from this follows B Q = J3.q 



3 Discussion and results 

This paper provides a rigorous framework for the description of linearized adiabatic la- 
grangian perturbations and stability of differentially rotating newtonian stars using semi- 
group theory. Problems of a previous framework by Dyson and Schutz are overcome and 
a basis for a rigorous analysis of the stability of such stars is provided. The spectrum of 
the oscillations is shown to coincide with the spectrum of an operator polynomial whose 
coefficients can be read off from the equation governing the oscillations about the equilib- 
rium configuration. Moreover, for the first time sufficient criteria for stability are given 
in form of inequalities for the coefficients of that polynomial. These show that a negative 
canonical energy of the star does not necessarily indicate instability. 

It is still unclear whether these criteria are strong enough to prove stability for realistic 
stars. On the other hand the second criterium has been sucessfully applied in the (on 
first sight seemingly unrelated case of the) stability discussion of the Kerr metric where 
the master equation governing perturbations is of the form ([]]), too.|| Another similarity 
of that case to the cases considered here is the fact that the corresponding operators C 
and B' 2 there are such that C — (1/A)B' 2 is positive whereas here this combination is 
semibounded as has been shown in DS. 

Also the determination of the spectrum of the operator polynomial C'—XB'+X 2 , A G C 
for some special case would be very useful. It is likely that this cannot be done for a 
physically relevant case. But it is also likely that the outcome to qualitative questions like 

• Does one have uniform stability in ml 

• Does a continuous part occur in the oscillation spectrum? 

only depends on structural properties of the operators C and B'. So from C probably 
only the highest order derivatives are relevant and details of the equation of state should 
be unimportant. From this point of view even the highly idealized case of a spherical 
background model with a truncated C along with a non constant velocity field v would 
be interesting to consider. 
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